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ABSTRACT 

Capitalizing on the observational advantage offered by its tiny M dwarf host, we present HST/WFC3 
grism measurements of the transmission spectrum of the super-Earth exoplanet GJ1214b. These are 
the first published WFC3 observations of a transiting exoplanet atmosphere. After correcting for a 
ramp-like instrumental systematic, we achieve nearly photon-limited precision in these observations, 
finding the transmission spectrum of GJ1214b to be flat between 1.1 and 1.7 /zm. Inconsistent with a 
cloud-free solar composition atmosphere at 8.2a, the measured achromatic transit depth most likely 
implies a large mean molecular weight for GJ1214b's outer envelope. A dense atmosphere rules out 
bulk compositions for GJ1214b that explain its large radius by the presence of a very low density 
gas layer surrounding the planet. High-altitude clouds can alternatively explain the flat transmission 
spectrum, but they would need to be optically thick up to 10 mbar or consist of particles with a range 
of sizes approaching 1 /im in diameter. 

Subject headings: planetary systems: individual (GJ 1214b) — eclipses — techniques: spectroscopic 



1. INTRODUCTION 
With a radius of 2.7 and a mass of 6.5 Mm, the 



transiting planet GJ1214b ( Charbonneau et al.||2009 ) is 
a member of the growing population of exoplanets whose 
masses and radii are known to be between those of Earth 
and Neptune ( sec Lcgcr ct ar'2009! ' Batalha et"aL||2011 
Lissauer et al. 1,2011; ^Winn et al.,|2(JTT r Among these 
exoplanets, most of which exhibit such shallow transits 
that they require ultra-precise space-based photometry 
simply to detect the existence of their transits, GJ1214b 
is unique. The diminutive 0.21 R© radius of its M dwarf 
stellar host means GJ1214b exhibits a large 1.4% tran- 
sit depth, and the system's proximity (13 pc) means the 
star is bright enough in the near infrared {H = 9.1) that 
follow-up observations to study the planet's atmosphere 
are currently feasible. In this work, we exploit this ob- 
servational advantage and present new measurements of 
the planet's atmosphere, which bear upon models for its 
interior composition and structure 



According to theoretica l studies (Seager et al. 



Rogers & Seager 



2010 
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Nettelmann et al. 2011 ), 



2007 



GJ1214b's 1.9 g cm bulk density is high enough to 
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require a larger ice or rock core fraction than the solar 
system ice giants but far too low to b e explained with an 
entirely Earth-like composition. Rogers & Seager (20101 
have proposed three general scenarios consistent with 
GJ1214b's large radius, where the planet could (i) have 
accreted and maintained a nebular H2 /He envelope atop 
an ice and rock core, (ii) consist of a rocky planet with 
an H2-rich envelope that formed by recent outgassing, or 
(iii) contain a large fraction of water in its interior sur- 
rounded by a dense H2-depleted, H20-rich atmosphere. 
Detailed the rmal evolution calculations by |Nettelmann| 



et al. (2011) disfavor this last model on the basis that 
it would require unreasonably large bulk water-to-rock 
ratios, arguing for at least a partial H2/He envelope, al- 
beit one that might be heavily enriched in H2O relative 
to the primordial nebula. 

By measuring GJ1214b's transmission spectrum, we 
can empirically constrain the mean molecular weight of 
the planet's atmosphere, thus distinguishing among these 
possibilities. When the planet passes in front of its 
host M dwarf, a small fraction of the star's light passes 
through the upper layers of the planet's atmosphere be- 
fore reaching us; the planet's transmission spectrum is 
then manifested in variations of the transit depth as a 
function of wavelength. The amplitude of the transit 
depth variations A_D(A) in the transmission spectrum 
scale as uh x 2HRp/R^, where tih is set by the opacities 
involved and can be 1-10 for strong absorption features, 
H is the atmospheric scale height, Rp is the planetary 
radiu s, and i^^, is the stellar radius (e.^ . Seager fc Sas- 

UOll. Because 



selov,i2000 Brown||2001 Hubbard et al. 



the scale height H is inversely proportional to the mean 
molecular weight /i of the atmosphere, the amplitude 
of features seen in the planet's transmission spectrum 
places strong constraints on the possible values of fi and, 
in parti cular, the hydrogen/heli um content of the atmo- 
sphere ( [Miller-Ricci et al.||2009[ ). 
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Indeed, detailed radiative transfer simulations of 



GJ1214b's atmosphere ( [Miller- Ricci fc Fortney] |2010 1 
show that a solar composition, H2-dominated atmo- 
sphere (/^ = 2.4) would show depth variations of roughly 
0.1% between 0.6 and 10 /im, while the features in an 
H20-dominated atmosphere (/i — 18) would be an order 
of magnitude smaller. While the latter of these is likely 
too small to detect directly with current instruments, the 
former is at a level that has regularly been measured with 



spectra of hot Jupiters (e.^ 


!;. Charbonneau et al. 


2002 


IPont et al.||2008 Sing et al. 


20ii 


). 




with 


Spectroscopic observations by 


Bean et al. 


(2010 



trum of GJ1214b to be featureless between 0.78-1.0 /xm, 
down to an amplitude that would rule out cloud-free H2- 
rich atmospheric models. Broadband Spizer Space Tele- 
scope photometric transit m easurements at 3.6 and 4.5 



fj,m by Desert et al. ( 2011a) sho wed a flat spectrum con- 
sistent wTtETBe^ret^r 1^010 1, as did high- resolution 



spectroscopy with MIK SPEC between 2.0 and 2.4 ^m by 



Crossfield et al. (2011). Intriguingly, the transit depth 



in A-band (2.2 lira) was measured from CFHT by Croll 
et al.| ( 20TT| to be 0.1% deeper than at other waveleiigths, 
which would imply a H2-rich atmosphere, in apparent 
contradiction to the other studies. 

These seemingly incongruous observations could po- 
tentially be brought into agreement if GJ1214b's at- 
mosphere were H 2-rich bu t significantly deplet ed in 
CH4 ( |Crossfield et ar||2011[ |Miller-Ricci Kcmpton et al. 



2011).^ In such a scenario, the molecular features that 



remain (predominantly H2O) would fit the CFHT mea- 
surement, but be unseen by the NIRSPEC and Spitzer 
observations. Explaining the flat VLT spectrum in this 
context would then require a broadband ha ze to smooth 



the spectrum at sho rter wavelengths (see Miller-Ricci 



Kempton et al. 2011). New observations by Bean et al 
( |201ip covering 0.6-0.85 fj,m and 2.0-2.3 /xm were again 
consistent with a flat spectrum, but they still could not 
directly speak to this possibility of a methane-depleted, 
H2-rich atmosphere with optically scattering hazes. 

Here, we present a new transmission spectrum of 
GJ1214b spanning 1.1 to 1.7 /im, using the infrared slit- 
less spectroscopy mode on the newly installed Wide Field 
Camera 3 (WFC3) aboard the Hubble Space Telescope 
(HST). Our WFC3 observations directly probe the pre- 
dicted strong 1.15 and 1.4 /im water absorption features 
in GJ1214b's atmosphere (Mil ler-Ricci fc Fortney||2010 1 
and provide a stringent constraint on the JI2 content of 
GJ1214b's atmosphere that is robust to non-equilibrium 
methane abundances and hence a definite test of the 
CH4-depleted hypothesis. The features probed by WFC3 
are the same features that define the J and H band win- 
dows in the telluric spectrum, and cannot be observed 
from the ground. 

Because this is the first published analysis of WFC3 ob- 
servations of a transiting exoplanet, we include a detailed 
discussion of the performance of WFC3 in this observa- 
tional regime and the systematic effects that are inherent 
to the instrument. Recent work on WFC3's predec essor 



NICMOS (Burke et al. 2010 



Gibson et al. 2011a) has 



highlighted the importance of characterizing instrumen- 
tal systematics when interpreting exoplanet results from 
HST observations. 
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Figure 1. A 512x100 pixel cutout of a typical WFC3 G141 grism 
exposure of the star GJ1214. The O"^ and l'^* order spectra are 
labeled, and the start of the 2"^* order spectrum is visible on the 
right. The location of the star in the direct images (not shown 
here) is marked with a circle. 

This paper is organized as follows: we describe our ob- 
servations in ^|2j our method for extracting spectropho- 
tometric light curves from them in ^|3j and our analy- 
sis of these light curves in ^|4 We present the resulting 
transmission spectrum and aiscuss its implications for 
GJ1214b's composition in f|5] and conclude in ^ 

2. OBSERVATIONS 

We observed three transits of GJ1214b on UT 2010 
October 8, 2011 March 28, and 2011 July 23 with the 
G141 grism on WFC3's infrared channel (HST Proposal 
#GO-12251, P.I. = Z. Berta), obtaining simultaneous 
multiwavelength spectrophotometry of each transit be- 
tween 1.1 and 1.7 /im. WFC3's IR channel consists of a 
1024 X 1024 pixel Teledyne HgCdTe detector with a 1.7 
/im cutoff that can be paired with any of 15 filters or 2 
low-resolution grisms (Dressel et al. 2010). Each expo- 
sure is compiled from multiple non-destructive readouts 
and can consist of either the full array or a concentric, 
smaller subarray. 

Each visit consisted of four 96 minute long HST orbits, 
each containing 45 minute gaps due to Earth occulta- 
tions. Instrumental overheads between the occultations 
are dominated by serial downloads of the WFC3 image 
buffer, during which all science images are transferred to 
the telescope's solid state recorder. This buffer can hold 
only two 16-readout, full-frame IR exposures before re- 
quiring a download, which takes 6 minutes. Exposures 
cannot be started nor stopped during a buffer download, 
so parallel buffer downloads are impossible for short ex- 
posures. 

Subject to these constraints and the possible readout 
sequences, we maximized the number of photons detected 
per orbit while avoiding saturation by gathering expo- 
sures using the 512 x 512 subarray with the RAPID 
NSAMP=7 readout sequence, for an effective integration 
time of 5.971 seconds per exposure. With this setting, 
four 12-exposure batches, separated by buffer downloads, 
were gathered per orbit resulting in an integration effi- 
ciency of 10%. Although the brightest pixel in the 1**' or- 
der spectrum reaches 78% of saturation during this expo- 
sure time, the WFC3's multiple non-destructive readouts 
enable the flux within each pixel to be estimated before 
the onset of significant near-saturation nonlinearities. 

A sub-region of a typical G141 grism image of GJ1214 
is shown in Fig. [T] The 512 x 512 subarray allows both 
the 0"^ and 1*^* order spectra to be recorded, and the 
2^ St Qj.(;[gj- spectrum to fall entirely within a single ampli- 
fier quadrant of the detector. The order spectrum 
spans 150 pixels with a dispersion in the x-direction of 
4.65 nm/pixel and a spatial full-width half maximum in 
the y-direction of 1.7 pixels (0.2"). The 0*^ order spec- 
trum is slightly dispersed by the grism's prism but is 
nearly a point source. Other stars are present in the 
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Figure 2. The mean out-of-transit extracted spectrum of GJ1214 
(black line) from all three HST visits, shown before (top) and af- 
ter (bottom) flux calibration. Individual extracted spectra from 
each visit are shown with their la uncertainties (color error bars). 
For comparison, the integrated flux from the PHOENIX model at- 
mosphere used to calculate the stellar limb darkening (see i ]4.2| is 
shown (gray lines) offset for clarity and binned to the WFCS pixel 
scale (gray circles). 

subarray's 68' x 61' field of view, but are too faint to 
provide useful diagnostics of systematic trends that may 
exist in the data. For wavelength calibration, we gath- 
ered direct images in the F130N narrow-band filter; the 
direct images' position relative to the grism images is 
also shown in Fig. [l] 

To avoid systematics from the detector flat-fie lds that 
have a quo ted precision no better than 0.5% (Pirzkal 
et al. 2011), the telescope was not dithered durmg any 
of the observations. We note that a technique called 
"spatial scanning" has been proposed to decrease the 
overheads for bright targets with WFC3, where the tele- 
scope nods during an exposure to smear the light along 
the cross-disp ersion direction, thus increasing the time 
to saturation ( [McCullough fc MacKe"nty||201li ). We did 
not use this mode of observation as it was not yet tested 
at the time our program was initiated. 



3. DATA REDUCTION 

The Python/PyRAF software package aXe was devel- 
oped to extract spectra from slitless grism obs ervations 
with WFC 3 and other Hubble instruments (Kiimmel 
et al.|2009 1, but it is optimized for extracting large num- 
bers of spectra from full frame dithered grism images. 
To produce relative spectrophotometric measurements of 
our single bright source, we opted to create our own ex- 
traction pipeline that prioritizes precision in the time 
domain. We outline the extraction procedure below. 

Through the extraction, we use calibrated 2- 
dimensional images, the "fit" outputs from WFC3's 
calwf 3 pipeline. For each exposure, calwf 3 performs the 
following steps: flag detector pixels with the appropriate 
data quality (DQ) warnings, estimate and remove bias 
drifts using the reference pixels, subtract dark current, 
determine count rates and identify cosmic rays by fitting 
a slope to the non-destructive reads, correct for photo- 
metric non-linearity (properly accounting for the signal 
accumulation before the initial "zeroth" read), and apply 
gain calibration. The resulting images are measured in 
e~s~^ and contain p er pixel uncertainty estim ates based 
on a detector model ( Kim Quijano et al.|2009 ). We note 



that calwf 3 does not apply fiat-field corrections when 
calibrating grism images; proper wavelength-dependent 
flat-fielding for slitless spectroscopy requires wavelength- 
calibrating individual sources and calwf 3 does not per- 
form this task. 

3.1. Interpolating over Cosmic Rays 

calwf 3 identifies cosmic rays that appear partway 
through an exposure by looking for deviations from a 
linear accumulation of charge among the non- destructive 
readouts, but it can not identify cosmic rays that appear 
between the zeroth and first readout. We supplement 
calwf 3's cosmic ray identifications by also flagging any 
pixel in an individual exposure that is > 6ct above the 
median of that pixel's value in all other exposures as a 
cosmic ray. Through all three visits (576 exposures), a 
total of 88 cosmic rays were identified within the extrac- 
tion box for the I''* order spectra. 

For each exposure, we spatially interpolate over cos- 
mic rays. Near the order spectrum, the pixel-to-pixel 
gradient of the point spread function (PSF) is typically 
much shallower along the dispersion direction than per- 
pendicular to it, so we use only horizontally adjacent 
pixels when interpolating to avoid errors in modeling the 
sharp cross-dispersion falloff. 

3.2. Identifying Continuously Bad Pixels 

We also mask any pixels that are identified as "bad 
detector pixels" (DQ=4), "unstable response" (DQ=32), 
"bad or uncertain fiat value" (DQ=512). We found that 
only these DQ fiags affected the photometry in a pixel 
by more than la. Other flags may have influenced the 
pixel photometry, but did so below the level of the pho- 
ton noise. In the second visit, we also identified one col- 
umn of the detector {x — 625 in physical pixel^]^ whose 
light curve exhibited a dramatically different systematic 
variation than did light curves from any of the other 
columns. This column was coincident with an unusually 
low-sensitivity feature in the flat-field, and we hypothe- 
size that the flat-field is more uncertain in this column 
than in neighboring columns. We masked all pixels in 
that column as bad. 

We opt not to interpolate over these continuously bad 
pixels. Because they remained flagged throughout the 
duration of each visit, we simply give these pixels zero 
weight when extracting ID spectra from the images. This 
allows us to keep track of the actual number of photons 
recorded in each exposure so we can better assess our 
predicted photometric uncertainties. 

3.3. Background Estimation 

In addition to the target, WFC3 also detects light from 
the diffuse sky background, which comes predominantly 
from zodiacal light and Earth-shine, and must be sub- 
tracted. We draw conservative masks around all sources 
that are visible in each visit's median image, inc luding 
GJ1214 and its electronics cross-talk artifact (see Viana 
fc Baggett||2010 ). We exclude these pixels, as well as all 



^ For ease of comparison with future WFC3 analyses, through- 
out this paper we quote all pixel positions in physical units as 
interpreted by SAOImage DS9, where the bottom loft pixel of a 
full-frame array would be (x,y) = (1,1). 
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Figure 3. Extracted properties of the 0**^ and l'^' order spectra as 
a function of time, including (a) the summed I''' order photometric 
hght curves; (b) the estimated sky background level; (c-g) the total 
flux, X and y position (measured relative to the reference pixel), 
and Gaussian widths in the x and y directions of the 0**^ order 
image; and (h-j) the y offset, cross-dispersion width, and slope of 
the I''' order spectrum. All three visits are shown and are denoted 
by the color of the symbols. The 45 minute gaps in each time series 
are due to Earth occultations, the 6 minute gaps are due to the 
WFC3 buffer downloads. 

pixels that have any DQ warning flagged. Then, to es- 
timate the sky background in e ach exposure, we scal e a 
master WFC3 grism sky image dKiimmel et"aL]|2011[ ) to 
match the remaining 70-80% of the pixels m each expo- 
sure and subtract it. We find typical background levels 
of 1 — 3 e~ s^^ pixel"^, that vary smoothly within or- 
bits and throughout visits as shown in Fig. p] (panel b). 
As a test, we also estimated the background level from 



a simple mean of the unmasked pixels; the results were 
unchanged. 

3.4. Inter-pixel Capacitance 

The normal calibration pipeline does not correct for 
the inter-pixel capacitance (IPC) effect, which effectively 
couples the fl ux recorded in adj acent pixels at about 
the 1% level ( |McCullough||2008[ ). We co rrect this ef- 
fect with a linear deconvolution algorithm (McCullough 



[2008'; 'Hilb ert fc McCullough||201l| ) , although we tind it 
makes little difference to the final results. 

3.5. Extracting the Zeroth Order Image 

The O*'' order image can act as a diagnostic for tracking 
changes in the telescope pointing and in the shape of the 
instrumental PSF. We select a 10 x 10 pixel box around 
the 0**^ order image and fit a 2D Gaussian to it with the 
X position, y position, size in the x direction, size in the y 
direction, and total flux allowed to vary (5 parameters). 

Time series of the 0**^ order x and y positions, sizes 
in both directions, and total flux are shown for all three 
visits in Fig. |3] (panels c-g) . Thanks to the dispersion by 
the grism's prism, the Gaussian is typically 20% wider in 
the X direction than in the y direction. Even though the 
throughput of the O"^ order image is a factor of 60 lower 
than the 1^* order spectrum, the transit of GJ1214b is 
readily apparent in the 0"^-order flux time series. 

3.6. Extracting the First Order Spectrum 

To extract the first order spectra, we first determine 
the position of GJ1214 in the direct image, which serves 
as a reference position for defining the trace and wave- 
length calibration of the I''' order spectrum. We adopt 
the mean position GJ1214 in all of the direct images as 
the reference position, which we measure using the same 
method as in extracting the 0*^ order image in {3.5 The 



measured {x,y) reference positions for the first, second, 
and third visits are (498.0, 527.5), (498.6, 531.1), and 
(498.9, 527.1) in physical pixels. 

Once the reference pixel for a visit is known, we use 
the coefflcien ts stored in the WFC3/ G141 aXe configu- 
ration filcpj ( Kuntschner et al. 2009), to determine the 
geometry of the 1"'' order trace and cut out a 30 pixel 
tall extraction box centered on the trace. Within this 
extraction box, we use the wavelength calibration coeffi- 
cients to determine the average wavelength of light that 
will be illuminating each pixel. We treat all pixels in the 
same column as having the same effective wavelength; 
given the spectrum's 0.5° tilt from to the x axis, errors 
i ntroduced by this simplifi cation are negligible. 



Kuntschner et al. (2008) used flat-fields taken through 
all narrow- band filters available on WFC3/IR to con- 
struct a fiat-field "cube" where each pixel contains 4 
polynomial coefficients that describe its sensitivity as a 
function of wavclcngtlf^ We use this fiat-field cube to 
construct a color-depcnHcnt flat based on our estimate of 
the effective wavelength illuminating each pixel, and di- 
vide each exposure by it. WFC3 wavelength calibration 
and flat-fielding is de scribed in detail in the aXe manual 
(IKiimmel et al. 120101). 



10 The aXe configuration file WFC3 . IR. G141 . V2 . . conf is avail- 
able through http://www.stsci.edu/hst/wfc3/ 

11 WFC3.IR.G141.flat.2.fits, through the same URL 



WFC3 Observations of GJ1214b 



5 



To calculate ID spectra from the flat-fielded images, 
we sum all the unmasked pixels within the extraction 
box over the j/-axis. To estimate the uncertainty in each 
spectral channel, we first construct a per-pixel uncer- 
tainty model that includes photon noise from the source 
and sky as well as 22 e^ of read noise, and sum these 
uncertainties, in quadrature, over the y-axis. We do not 
use the calwf 3-estimated uncertainties; they include a 
term propagated from the uncertainty in the nonlinearity 
correction that, while appropriate for absolute photome- 
try, would not be appropriate for relative photometry. In 
each exposure, there are typically 1.2 x 10^ e~ per single- 
pixel spectral channel and a total of 1.5 x 10^ e~ in the 
entire spectrum. Fig. Is] (panel a) shows the extracted 
spectra summed over all wavelengths as a function of 
time, the "white" light curve. 

For diagnostics' sake, we also measure the geometri- 
cal properties of the I''* order spectra in each exposure. 
We fit ID Gaussians the cross-dispersion profile in each 
column of the spectrum and take the median Gaussian 
width among all the columns as a measurement of the 
PSF's width. We fit a line to the location of the Gaus- 
sian peaks in all the columns, taking the intercept and 
the slope of that line as an estimate of the y-offset and 
tilt of the spectrum on the detector. Time series of these 
parameters are shown in Fig. |3] (panels h-j). 

3.7. Flux Calibration 

For the sake of display purposes only (see Fig. [2]) , we 
flux calibrate each visit's median, extracted, ID spec- 
trum. Here we have interpolated over all bad pixels 
within each visit (contrary to the discussion in the ^3.2), 
and plotted the weighted mean over all three visits. The 
calibration uncertainty for the G141 sensitivity curve 
(iKuntschner et al.||2011|) is quoted to be 1%. 



3.8. Times of Observations 

For each exposure, we extract the EXPSTART keyword 
from the science header, which is the Modified Julian 
Date at the start of the exposure. We correct this to 
the mid-exposure time using the EXPTIME keyword, and 
convert it to the Barycentric Julian Date in the Barycen- 
tri c Dynamical Time sta ndard using the code provided 
by [Eastman et al.|p010l). 



4. ANALYSIS 

In this section we descr ibe our method for estimating 
parameter uncertainties (i ]4.1[ ) and our s trate gy for mod- 
eling GJ1214's stellar limE~darkening (jj42|). Then, af- 
ter ident ifyin g the dominant systematics m WFC3 light 
curves (^4.3) and describing a method to correct them 
(j|4.4[), we present our fi ts to the light curves, both 



summed over wavelength (j |4.5[ ) and spectroscopically re- 
solved (^4.6|. We also present a fruitless sear ch for tran- 
siting satellite companions to GJ1214b (^4.7). 



4.1. Estimating Parameter Distributions 

Throughout our analysis, we fit different WFC3 light 
curves with models that have different sets of parame- 
ters, and draw conclusions from the inferred probability 
distributions of those parameters; this section describes 
our method for characterizing the posterior probability 
distribution for a set of parameters within a given model. 



We use a Markov Chain Monte Carlo (MCMC) method 
with the Metropolis-Hastings algorithm to explore the 
posterior probability density function (PDF) of the 
model parameters. This Bayesian technique allows us 
to sample from (and thus infer the shape of) the prob- 
ability distribution of a model's parameters given both 
our data and our prior knowledge about the parameters 
(for r eviews, see Ford |2005[ Gregory 2005 Hogg et al. 
20101. Briefly, the algorithrn starts a chain with an ini- 
tial set of parameters (Mj=o) and generates a trial set of 
parameters (M'j+i) by perturbing the previous set. The 
ratio of posterior probability between the two parameter 
sets, given the data D, is then calculated as 



P(M^,+i|D) ^ P(D|M',+i) PiM',+,) 
F(M,|D) P(D|M,) P(M,) 



(1) 



where the first term (the "likelihood") accounts for the 
information that our data provide about the parameters 
and the second term (the "prior" ) specifies our externally 
conceived knowledge about the parameters. If a random 
number drawn from a uniform distribution between 
and 1 is less than this probability ratio, then Mj+i is set 
to M'j-|_i; if not, then Mj+i reverts to Mj. The process 
is iterated until j is large, and the resulting chain of 
parameter sets is a fair sample from the posterior PDF 
and can be used to estimate confidence intervals for each 
parameter. 

To calculate the likelihood term in Eq. [T] we assume 
that each of the N fiux values di is drawn from a uncorre- 
lated Gaussian distribution centered on the model value 
mi with a standard deviation of stXi, where ai is the the- 
oretical uncertainty for the flux measurement based on 
the detector model and photon statistics and s is a pho- 
tometric uncertainty rescaling parameter. Calculation of 
the ratio in Eq. [T] is best done in logarithmic space for 
numerical stability, so we write the likelihood as 



1 



lnP(D|M) = -A^lns - —^x + constant (2) 



where 



2s2 



N , , X 2 

d, - m, 



(3) 



and we have only explicitly displayed terms that depend 
on the model parameters. Including s as a model pa- 
rameter is akin to rescaling the uncertainties by exter- 
nally modifying ai to achieve a reduced of unity, but 
enables the MCMC to fit for and marginalize over this 
rescaling automatically. Unless otherwise stated for spe- 
cific parameters, we use non-informative (uniform) priors 
for the second term in Eq. [T] We use a Jeffreys prior on s 
(uniform in In s) which is the least informative, although 
the results are practically indistinguishable from prior 
uniform in s. 

When g enerating each new tri al parameter set M'j+i, 
we follow Dunkley et al. (2005) and perturb every pa- 
rameter at once, drawing the parameter jumps from a 
multivariate Gaussian with a covariance matrix that ap- 
proximates that of the parameter distribution. Doing so 
allows the MCMC to move easily along the dominant lin- 
ear correlations in parameter PDF, and greatly increases 
the efficiency of the algorithm. While this procedure may 
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Table 1 

White Light Curves from WFC3/G141 



Timc^ 


Relative Flux° 


Uncertainty 


Sky 




Qth_Yc 








Is'-B^ 


Ist-Slope'^ 


Visit 


(BJDtdb) 






{e-/s) 


(pix) 


(pix) 


(pix) 


(pix) 


(pix) 


(pix) 


(pix/pix) 




2455478.439980 


0.99381 


0.00031 


1.9546 


-187.830 


-0.474 


0.790 


0.613 


-0.075 


0.7484 


0.00921 


1 


2455478.440270 


0.99713 


0.00031 


1.9938 


-187.839 


-0.481 


0.792 


0.616 


-0.082 


0.7490 


0.00925 


1 


2455478.440559 


0.99787 


0.00031 


1.9718 


-187.846 


-0.484 


0.789 


0.615 


-0.076 


0.7486 


0.00914 


1 


2455478.440848 


0.99958 


0.00032 


1.8808 


-187.827 


-0.490 


0.792 


0.614 


-0.087 


0.7476 


0.00917 


1 


2455478.441138 


0.99989 


0.00032 


1.9379 


-187.844 


-0.490 


0.785 


0.612 


-0.085 


0.7492 


0.00921 


1 



Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is 
shown here for guidance concerning its form and content. 
^ Mid-exposure time. 

^ Normalized to the median flux level of the out-of-transit observations in each visit. 
Position measured relative to the Gaussian center of each visit's direct image. 

Gaussian width of the 0'^ or 1^* order spectra in the horizontal (A) or vertical (B) direction. 

" Slope of the 1^' order spectrum. 



seem circular (if we knew the covariance matrix of the pa- 
rameter distribution, why would we need to perform the 
MCMC?), the covariance matrix we use to generate trial 
parameters could be a very rough approximation to the 
true shape of the parameter PDF but still dramatically 
decrease the computation time necessary for the MCMC. 

To obtain an initial guess for p arameters (Mj^p) , 
we use the MPFIT implementation (Markwardt 20091 
of the Levenberg-Marquardt (LM) method to maximize 
lnP(M|D). This would be identical to minimizing 
in the case of flat priors, but it can also include con- 
straints from more informative priors. The LM fit also 
provides an estimate of the covariance matrix of the pa- 
rameters, which is a linearization of the probability space 
near the best-fit. We use this covariance matrix esti- 
mate for generating trial parameters in the MCMC, and 
with it, achieve parameter acceptance rates of 10-40% 
throughout the following sections. As expected, when 
fitting models with flat priors and linear or nearly-linear 
parameters (where the PDF should well-described by 
a multivariate Gaussian), the LM covariance matrix is 
iden tical to tha t ultimately obtained from the MCMC 
(see Sivia||1996| for further discussion). 

MCMC chains are run until they contain 1.25 x 10^ 
points. The first 1/5 of the points are ignored as "burn- 
in" , leaving 1 x 10^ for parameter estimation. Correlation 
lengths for the parameters in the MCMC chains are indi- 
cated throughout the text; they are typically of order 10 
points. A chain with such a correlation length effectively 
contains 1 x 10^/10 = 1 x lO'' independent realizations 
of the posterior PDF. We quote confidence intervals that 
exclude the upper and lower 16% of the marginalized dis- 
tribution for each parameter (i.e. the parameter's central 
68% confidence interval), using all 1x10^ points in each 
chain. 

4.2. Modeling Stellar Limb Darkening 

Accurate modeling of the WFC3 integrated and spec- 
troscopic transit light curves requires careful consid- 
eration of the stellar limb-darkening (LD) behavior. 
GJ1214b's M4.5V stellar host is so cool that it exhibits 
weak absorption features due to molecular H2O. Because 
inferences of the planet's apparent radii from transit light 
curves depend strongly on the star's limb-darkening. 



which is clearly influenced by H2O as an opacity source, 
inaccurate treatment of limb- darkening could potentially 
introduce spurious H2O features into the transmission 
spectrum. 

If they were sufficiently precise, transit light curves 
alone could simultaneously constrain both the star's mul- 
t i wavelengt h limb- darkenin g behavior and the pl anet ' s 
multiwavelength radii (e.g. Knutson et al. 2007b). For 
less precise light curves, it is common practice to fix the 
limb-darkening to a theoretically calculated law, even if 
this may undere stimate the uncertainty in the pla netary 
parameters (see Burke et al. 2007 Southworth 2008j). 
Given the quality of our data, we adopt an mterniediate 
solution where we allow the limb-darkening parameters 
to vary in our fits, but with a Gaussian prior centered on 
the theoretical values (e.g. Bean et aL]|2010 ). 

We model the star GJl214's limb-darkening behav- 
ior with a spherically symmetric PHOENIX atmosphere 
(Hauschildt et al. 19991, assuming stellar parameters of 
I'es = 3026 A, logq^5, and [M/H] = ( [Charbonneau 
et al. 2009|). As shown in Fig. [2] the integrated flux 



from the PHOENIX model is in good qualitative agree- 
ment with the low-resolution, calibrated WFC3 stellar 
spectrum of GJ1214. From this model, we calculate 
photon-weighted average intensity profiles for the inte- 
grated spectrum and for each of the individual wave- 
length bins, using the WFC3 grism sensitivity curve and 
the PHOENIX model to estimate the photon counts. In 
the spherical geometry of the PHOENIX atmospheres 
the characterization of the actual limb (defined as /i = 0, 
see below) is not straightforward, as the model extends 
beyond the photosphere into the optically thin outer at- 
mosphere. The result is an approximately exponentially 
declining intensity profile fr om the outermost layers, that 
I Claret &: Hauschildt (2003) found not to be easily repro- 
duced by standard iini5~darkening laws for plane-parallel 
atmospheres. These authors suggest the use of "quasi- 
spherical" models by ignoring the outer region. In an 
extension of this concept, we set the outer surface of the 
star to be where the intensity drops to e~^ of the central 
intensity, and measure /i — cos 9 (where 6 is the emission 
angle relative to the line of sight) relative to that outer 
radius. 

We derive coefficients for a square-root limb-darkening 



WFC3 Observations of GJ1214b 



7 



V 1.02 
§ 1.01 
S 1.00 
I 0.99 
|0.98 

V1.02 
§ 1.01 
f 1.00 
§ 0.99 
|0.98 

C 1.02 
§ 1.01 
f 1.00 
§ 0.99 
|0.98 

V 1.02 

§ 1.00 
§ 0.99 
|o.98 

V 1.02 

§ 1-01 
f 1.00 
§0.99 
io.98 



-i---f----t- 



4 



J ^ 



% ^> i- 

til 




12 3 4 

Time from Start of Each 12-exDosure Botch ('minutes) 

Figure 4. Single-pixel light curves within each 12-exposure batch 
following a buffer download {gray points), shown for different mean 
pixel illuminations. Pixel light curves have been normalized to the 
first exposure within each batch, and plotted with small random 
horizontal offsets for clarity. Only data from the first HST visit, 
which exhibited the smallest pointing drifts (see Fig.jsJ, are shown. 
Error bars show the mean and its standard error for each time point 
and each illumination. An exponential ramp begins is present in 
pixels with a mean recorded fluence greater than 30,000 and 40,000 
e~ (50% of the detector full well). Note, the nominal fluences 
quoted here do not include charge accumulated during detector 
flushing and initial readout (see text). 

law for each of these average intensity profiles using least- 
squares fitting. In this law, the intensity relative to the 
center of the star is given by 



i-c{i- ^l)-d(^-^), 



(4) 



where c and d are the two coefficients of the fit. We 
chose a square-root law over the popular quadratic law 
because it gave noticeably better approximations to the 
PHOENIX intensity profiles, while still having few enough 
free parameters t hat they can be par tially inferred from 
the data. Indeed, van Hamme ( 1993 1 found the square- 
root law to be generally preferable to other 2-parameter 
limb-darkening laws for late-type stars in the near-IR. 
The square-root law matches the theoretical intensity 
profile nearly as well as the full nonlinear 4-parameter 
law introduced by Claret (2000) for the models we use 
here. 

4.3. Light Curve Systematics 

The summed light curve shown in Fig. [s] (panel a) ex- 
hibits non-astrophysical systematic trends. The most ob- 
vious of these are the sharply rising but quickly saturat- 
ing "ramp" -like features within each batch of 12 expo- 



sures between buffer downloads. To the eye, the ramps 
are very repeatable; the flux at the end of all batches 
asymptotes to nearly the same level. The amplitude of 
the ramp is 0.4% from start to finish for most batches, 
except for the first batch of each orbit, where the ramp 
is somewhat less pronounced. 
These ramps are reminiscent of those seen in high- 



cadenc e Spitzer light curves at 8 and 16 /xm (e.g. Dem- 



ing et a l. 2006; Kn utson et al.|2007a Charbonneau et aT 



2008[ ) which |Agol et al. ( |2010p recently proposed may 
be due to "charge trapping" within the detector pixels. 
In their toy model, charge traps within each pixel be- 
come filled throughout an exposure and later release the 
trapped charge on a finite timescale, thereby increasing 
the pixel's dark current in subsequent exposures. The 
model leads to exponential ramps when observing bright 
sources as the excess dark current increases sharply at 
first but slows its increase as the population of charge 
traps begins to approach steady state. We note this 
model also leads to after-images following strong expo- 
sures, i.e., persistence. 

WFC3 has been known since its initial ground-testing 
to exhibit strong persistence be h avior ([McCul lough &:[ 

I 



DeustuaJ [20081 
have proposec 



Long et al. 2010). Smith et al. (20083 



that persistence m 1.7 /zm cutoti" HgCdle 



detectors like WFC3 is likely related to charge trapping. 
Measurements (McCullough & Deustua 2008) indicate 
that WFC3's persistence may be of the right order of 
magnitude (on < 1 minute timescales) to supply the 
roughly 50 e" s~^ pixel" "'^ in the brightest pixels that 
would be necessary to explain the observed several mil- 
limagnitude ramp, although persistence levels and de- 
cay timescales can depend in com plicated ways on th e 
strength of previous exposures (see Smith et al.||20 08"b) . 

We were aware of this persistence issue before our ob- 
servations and made an effort to control its effect on our 
light curves. When we planned the timing of the expo- 
sures, we attempted to make the illumination history of 
each pixel as consistent as possible from batch to batch 
and orbit to orbit. In practice, this means we gathered 
more direct images than necessary for wavelength cali- 
bration to delay some of the grism exposures. 

Whether or not the ramps are caused by the charge 
trapping mechanism, they are definitely dependent on 
the illumination that a pixel receives. To demonstrate 
this, we construct light curves for each individual pixel 
over the duration of every out-of-transit 12-exposure 
batch that follows a buffer download and normalize each 
of these pixel light curves to the first exposure in the 
batch. Fig. H] shows the normalized pixel light curves, 
grouped by their mean recorded fluence. Because it takes 
a finite time to read the subarray (0.8 seconds) and reset 
the full array (2.9 seconds), we note that each exposure 
actually collects 60% more electrons than indicated by 



these nominal, recorded fluences (see Long et al. 2011) 



The appearance of the ramp clearly becomes more pro- 
nounced for pixels that are more strongly exposed. 

Buried beneath the ramp features, the summed light 
curve exhibits subtler trends that appear mostly as orbit- 
long or visit-long slopes with a peak-to-peak variation of 
about 0.05%. These are perhaps caused by slow drifts 
in pointing and focus (telescope "breathing" ) interacting 
with sensitivity variations across the detector that are 
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not perfectly corrected by the flat field. 

4.4. Correcting for Systematics 

Fortunately, these systematics are extremely repeat- 
able between orbits within a visit; we harness this fact 
when correcting for them. We divide the in-transit 
orbit of any photometric timeseries, either the white 
light curve or one of the spectroscopically resolved light 
curves, by a systematics correction template constructed 
from the two good out-of-transit orbits. This template is 
simply the weighted average of the fluxes in the out-of- 
transit orbits, evaluated at each exposure within an orbit. 
It encodes both variations in the effective sensitivity of 
the detector within an orbit and the mean out-of-transit 
flux level. 

When performing the division, we propagate the tem- 
plate uncertainty into the photometric uncertainty for 
each exposure, which typically increases it by a factor 
of y/1 + 1/2 = 1.22. This factor, ahhough it may seem 
like an undesired degradation of the photometric preci- 
sion, would inevitably propagate into measurements of 
the transit depth whether we performed this correction 
or not, since Rp/Ri, is always measured relative to the 
out-of-transit flux, which must at some point be inferred 
from the data. 

Throughout this work, we refer to this pro- 
cess of dividing by the out-of-transit orbits as the 
divide-oot method. Because each point in the sin- 
gle in-transit orbit is equally spaced in time between 
the two out-of-transit exposures being used to correct 
it, the divide-oot method also naturally removes the 
0.05% visit- long slope seen in the raw photometry. As 
we show in |4.5( when applied to the white light curves, 
the divide-oot treatment produces uncorrelated Gaus- 
sian residuals that have a scatter consistent with the pre- 
dicted photon uncertainties. 

Unlike decorrelation techniques that have often been 
used to correct systematics in HST light curves, the 
divide-oot method does not require knowing the re- 
lationship between measured photometry and the phys- 
ical state of the camera. It does, however, strictly re- 
quire the systematics to repeat over multiple orbits. The 
divide-oot method would not work if the changes in 
the position, shape, and rotational angle of the f^'-order 
spectrum were not repeated in the other orbits in a visit 
or if the cadence of the illumination were not nearly iden- 
tical across orbits. In such cases, the Ga ussian process 
method proposed by Gibson et al. (2011b) may be a use- 
ful alternative, and one that would appropriately account 
for the uncertainty involved in the systematics correc- 
tion. 



4.5. White Light Curve Fits 

Although the main scientific result of this paper is de- 
rived from the spectroscopic light curves presented in 
^4.6 we also analyze the light curve summed over all 
wavelengths between 1.1 and 1.7 /im. We use these 
white light curves to confirm the general system prop- 
erties found in previous studies and quantitatively inves- 
tigate the instrumental systematics. 

We fit an analytic, limb-darke ned transit light curve 
model ( [Mandel fc Agol| |2002[ ) to the divide-oot- 
corrected white light curves. Only the in-transit orbits 



were fit; after the divide-oot correction, the two out- 
of-transit orbits contain no further information. Also, 
because the in-transit orbit's flux has already been nor- 
malized, we fix the out-of-transit flux level to unity 
in all the fits. Throughout, we fix the planet's pe- 
riod to P = 1.58040481 days a nd mid-transit tim e to 
Tc = 2454966.525123 BJDtdb ( |Bean et al.||2011[ ), the 
orbital ec centricity to e = 0. and the stellar mass to 
O.157M0 ( [Charbonneau et al.||2009"l ). 



4.5.1. Comhmed White Light Curve 

First, we combine the three visits into a single light 
curve, as shown in Fig. [5j and fit for the following pa- 
rameters: the planet-to-star radius ratio (i?p/i?*), the 
total transit duration between first and fourth contact 
(ti4), the stellar radius (i?*), and the two coefficients c 
and d of the square-root limb-darkening lavij^ Previous 
studies have found no significant t ransit timing variations 
for the G J 1214b system (Charbonneau et al.||2009; Sada 



eral] |20l0'; 'Be an et al |[2010( i Carier ei al.||mi| ipesert 



eF^]]2011a; Kundurthy et al ||2011[ |Berta et al.||2011[ 
CroU et ar]p01l[ ), so we fix the time of mid-transit for 
each vis itto be that predict ed by the linear ephemeris. 

As in Burke et al. (2007), we use the parameters tu 
and i?+ to ensure quick convergence of the MCMC be- 
cause correlations among these parameters are more lin- 
ear than for the commonly fit impact parameter (6) and 
scaled semi-major axis {a/ R^). Because nonlinear trans- 
formations between parameter pairs will deform the hy- 
pervolume of parameter space, we include a Jacobian 
term in the priors in Eq. fl] to ensure unifo rm priors for 
the physical param eters Rp, R^,, and i (see Burke et al 



2007 Carter et al. 2008, for detailed discussions]. Fot 



the combined light curves, the influence of this term is 
practically negligible, but we include it for completeness. 
In the MCMC chains described in this section, all pa- 
rameters have correlation lengths of 6-13 points. 

Initially, we perform the fit with limb-darkening coeffi- 
cients c and d without any priors from the PHOENIX atmo- 
sphere model, enforcing only that < c -I- d < 1, which 
ensures that the star is brighter at its center (fi = 1) 
than at its limb = 0). Interestingly, the quantity 
(c/3 + d/5), which sets the integral of /(^) over the 
stellar surface, defines the line along which c and d are 



Irwin et al.| 


201l|) 


commonly quoted 



For quadratic limb- darkening, the 



2ui + U2 combination ( Holman et al 
20061 has the same physical meaning, ihe integral of 
1 {ji) can be thought of as the increase in the central 
transit depth over that for a constant-intensity stellar 
disk, so it makes sense that it is well-constrained for 
nearly equatorial transiting systems like GJ1214b. Plan- 
ets with higher impact parameters do not sample the full 
range of < /j, < 1 during transit, leading to correspond- 
ingly weaker limb-darkening const raints that can be de - 
rived from their light curves (see Knutson et al.[|2011 ). 
We quote confidence intervals for the linear combination 
(c/3 -I- d/5) and one orthogonal to it in Table [2j along 
with rest of the parameters. 

Heartened by finding that when they are allowed to 
vary freely, our inferred white-light limb-darkening coef- 

The square-root law is a special case of the 4-parameter law 
and straightforward to include in the ^Mandel &: Agol_( 2002 ) model. 
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after the correction has been 



model-ramp 



model-ramp 
Visit 1 
Visit 2 
Visit 3 
combined 



Figure 5. The white light curve of GJ1214b's transits before {top panels) and after (middle pane ls) re mov ing th e instrumental systematics 
using the divide-oot (left) and model-rcunp (right, with offsets for clarity) methods describ ed in jjX4|and j]4.5.4| A transit model that was 
fit to the divide-oot-corrected light curve, constrained to the values of a/i?* and b used by |Bean et al.| | [20iu[ l, is shown (gray lines), along 
with residuals from this model (bottom panels). In the left panels, the out-of-transit orbits are not shown ai 
applied, because they contain no further information. 

ficients agree to la to those derived using the PHOENIX 
stehar model, we perform a second fit that includes the 
PHOENIX models as informative priors. For this prior, 
we say P(M) in Eq. [1] is proportional to a Gaussian 
with (c/3 + d/5) = 0.TI892 ± 0.018 and (c/5 - d/3) = 
-0.431 ±0.032, which is centered on the PHOENIX model. 
To set the la widths of these priors, we start by vary- 
ing the effective temperature of the star in the PHOENIX 
model by its 130K uncertainty in either direction, and 
then double the width of the prior beyond this, to ac- 
count for potential systematic uncertainties in the atmo- 
sphere model. The results from the fit with these LD 
priors are shown in Table [2j 

The photometric noise rescaling parameter s is within 
10% of unity, implying that the 376 ppm achieved scat- 
ter in the combined white light curve can be quite well- 
explained from the known sources of uncertainty in the 
measurements, predominantly photon noise from the 
star. As shown in Fig. [6] for the divide-oot-corrected 
light curves, the autocorrelation function (ACF) of the 
residuals shows no evidence for time-correlated noise. 
Likewise, the scatter in binned divide-oot residuals de- 
creases as the square-root of the number of points in a 
bin, as expected for uncorrelated Gaussian noise. If there 
are uncorrected systematic effects remaining in the data, 
they are below the level of the photon noise over the 
time-scales of interest here. 




5 10 15 5 10 15 

Log [top] or Bin Size [bottom] (dotopoints) 

Figure 6. For the transit model in Fig. [s] and both types of sys- 
tematics treatments, the autocorrelation function of the residuals 
(ACF; top) and the scatter in binned residuals as a function of 
bin size (bottom). The residuals from the combined light curve 
are shown (black points), as well as the individual visits (colorful 
points). The expectations from uncorrelated Gaussian noise (0 in 
the top, <x 1/\/N in the bottom) are overplotted (dashed lines). 



durthy et al. 2011 
sured A'J 



Berta et al. 
s show no evidence : 



4.5.2. Individual White Light Curves 

To test for possible differences among our WFC3 vis- 
its, we fit each of the three divide-oot-corrected white 
light curves individually. In addition to i?p/i?^, tu, and 
Ri,, we also allow ATc (the deviation of each visit's mid- 
transit time from the linear ephemeris) to vary freely. 
We allow c and d to vary, b ut enfo rce the same PHOENIX- 
derived priors described in §4.5.1 



Table [3] shows the results, which are c onsistent with 
each other and wit h other o bservations fCharbonneau' 
eTaL] [20091 |Bean et al., 20101 [Carter et al.,2011, ,Kun-| 



2011| ). The three mea- 
or transit timing varia- 
tions. The uncertainties for the parameters ti4, R^, and 
ATc are noticeably largest in the first visit; this is most 
likely because the first visit does not directly measure 
the timing of either I*'* or 4*^ contact, on which these 
parameters strongly depend. Additionally, whereas the 
correlation lengths in the MCMC chains for these pa- 
rameters in the two visits that do measure l*''/4*^ con- 
tact and for i?p/i?^, c, and d in all three visits are small 
(10-30 points), the correlation lengths for ti4, i?^, and 
ATc in the first visit are very large (300-400 points), in- 
dicating these weakly constrained parameters are poorly 
approximated by the MPFIT-derived covariance matrix. 
On account of the large correlation lengths for these pa- 
rameters, we ran the JVICMC for the first visit with a 
factor of 10 more points. In each of the three visits, the 
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Table 2 

Transit Parameters Inferred from the Combined White Light 
Curve 



Parameter 



No LD Prior 



LD Prior^ 



ii4 (days) 

(Rq) 
a/-R* 
i (°) 
b 

c/3 + d/5 
c/5 - d/3 

predicted RMS 
achieved RMS 



0.03624 ± 0.00013 

U.ZUlt_Q QQ25 

10.0U_Q 29 

89.3 ±0.4 

o.isin? 

U.liOO_Q QQQg 

0.096 ± 0.008 

-,+0.22 
-0.14 



0.03620 ± 0.00012 
201+'^ °°'' 

10.01_Q 29 

89.3+«-^ 



0.19 



-0.3 
-0.08 



-0.52^ 



337 ppm 
373 ppm 
1.12 ±0.07 



-0.11 

0.1160 ±0.0005 
0.095 ±0.007 
-0.433 ±0.032 

337 ppm 
376 ppm 
1.12 ±0.07 



^ The Gaussian limb-darkening priors of (c/3 + d/5) = 
0.0892 ± 0.018 and c/5 - d/3 = -0.4306 ± 0.032 were 
derived from PHOENIX stellar atmospheres, as described 
in the text. 

^ Confidence intervals on Rp/R^, do not include the 
ARp/Ri, = 0.0006 systematic uncertainty due to stellar 
variability (see text). 



Table 3 

Transit Parameters Inferred from Individual White Light Curves 



Parameter 



Visit 1 



Visit 2 



Visit 3 



ATc (days)'' 
ti4 (days) 

R* (R©) 
a/Ri, 

i n 

b 

Rp/R* 

predicted RMS 
acheived RMS 



-O.OOOllo:oo20 0.00028 ± 0.00031 0.0002 ± 0.0004 



0.037 
0.211 
14 

88.9 ±0.7 
0.27 ± 0.17 
0.1164 



+0.004 
0.003 
+0.021 
0.014 
-1 
-1 



+0.0009 
0.0008 



0.0357 ± 0.0007 
0.200j 
15.4:! 
89.2 ± 0.6 
0.21 ± 0.15 
0.1159^ 



f0.008 
-0.006 
.1+0.5 
-0.6 



-j+O.OOll 
-0.0009 



0.0369 ±0.0010 

U.^14_0.011 
14.4 ± 0.9 
88. 5^ 
0.381 
0.1175^ 



^+0.8 
-0.7 
s+0.13 
-0.20 
-+0.0011 
-0.0012 



337 ppm 
343 ppm 

^■'^'-o.ii 



337 ppm 
360 ppm 

1 1 9+0.13 
^■^^-0.11 



337 ppm 

366 ppm 

1 14+0" 
^•^^-0.11 



^ Offset between the observed mid-transit time and that cal- 
culated from the linear ephemeris with P = 1.58040481 and 
Tc=2454966.525123 BJDtdb- 

uncertainty rescaling parameter s is slightly above but 
consistent with unity, indicating the photometric scatter 
is quite well explained by known sources of noise. 



4.5.3. Stellar Vanabtlity 



GJ1214 is known to be variable on 50-100 day 
timescales with an ampl itude of 1% in the MEarth band- 
Charbonneau et al.||2009] [Berta 



pass (715-1000 nm; see 



et al. 12011 1 . To gauge the impact of stellar variability in 



the wavelengths studied here, we plot in Fig. [7] the rela- 
tive out-of-transit flux as measured by our WFC3 data. 
For each HST visit, we have three independent measure- 
ments of this quantity: the F130N narrow-band direct 
image, the O^'^-order spectrum, and the I'^'-order spec- 
trum. Consistent variability over these measurements 
that sample different regions of the detector within each 
visit would be difficult to reproduce by instrumental ef- 
fects, such as flat-fielding errors. In Fig. [Tj GJ1214 ap- 
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Figure 7. The relative out-of-transit (O.o.T.) flux for each HST 
visit, measured independently from three different groups of im- 
ages; the summed I'^'-order spectrum, the O'^'-order image, and 
the narrow-band direct image, each normalized to its mean. Error 
bars denote the standard deviation of the out-of-transit measure- 
ments within each visit; they do not include the 0.5% uncertainty 
in the detector flat-field. The narrow-band measurements sample 
fewer photons, thus their larger uncertainties. 

pears brighter in the first visit than in the last two visits, 
with an overall amplitude of variation of about 1%. 

This 1% variability, if caused by unoccultcd spots 
on the stellar surface, should lead to variations in 
the inferred planet -to-star radius ra tio on the order of 
ARp/R^, = 0.0006 ( [Berta et al.|2011[ ). This is larger than 
the formal error on Rp/ Hi, troni the combined white light 
curve (Tab. [2]), and must be considered as an important 
systematic noise floor in the measurement of the abso- 
lute, white-light transit depth. We do not detect this 
variability in the individually measured transit depths 
(Tab. pi) because it is smaller than the uncertainty on 
each. Most importantly, while the spot-induced variabil- 
ity influences the absolute depth at each epoch, its effect 
on the relative transit depth among wavelengths will be 
much smaller and not substantially bias our transmission 
spectrum estimate. 

4.5.4. Modeling Instrumental Systematics 

Before calculating GJ1214b's transmission spectrum, 
we detour slightly to use the white light curves's high 
photometric precision to investigate the characteristics 
of WFC3's instrumental systematics. Rather than cor- 
recting for the instrumental systematics with the simple 
non-parametric divide-oot method, in this section we 
describe them with an analytic model whose parameters 
illuminate the physical processes at play. We refer to this 
treatment as the model-rainp method. 

In this model, we treat the systematics as consisting 
of an exponential ramp, an orbit-long slope, and a visit- 
long slope. We relate the observed flux (Fobs) to the 
systematics-free flux (Fcor) by 



obs 



{C + ytvis + SWb) (1 - i?e-(*-'-^-)/^) (5) 



where ivis is time within a visit (= at the middle of each 
visit), torb is time within an orbit {— at the middle of 
each orbit), that is time within a batch (= at the start 
of each batch), t is a ramp timescale, and the term 



D 




for the 1^' batch of an orbit 
for the other batches 



(6) 



allows the exponential ramp to be delayed slightly for the 
first batch of an orbit. 

The ex ponential form aris es out of the toy model pro- 
posed by Agol et al. (20101, where a certain volume of 
the detector pixels lias the ability to temporarily trap 
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charge carriers and later release them as excess dark cur- 
rent. In quick series of sufficiently strong exposures, the 
population of charge traps approaches steady state, cor- 
responding to the flattening of the exponential. Judging 
by the appearance of the ramp in the 2"'*-4*'' batches of 
each orbit, the release timescale seems to be short enough 
that the trap population completely resets to the same 
baseline level after each 6 minute buffer download (dur- 
ing which the detector was being continually flushed each 
2.9 seconds). Compared to these batches, the batch 
of each orbit appears to exhibit a ramp that is either 
weaker, or as we have parameterized it with the Df, term, 
delayed. We do not explain this, but we hypothesize that 
it relates to rapid changes in the physical state of the 
detector coming out of Earth occultation affecting the 
pixels' equilibrium charge trap populations. 

The visit-long and orbit-long slopes are purely descrip 
tive terms (as in Brown et al. 2001 



Nutzman et al.|2(ill[ ), but relate to p hysical processes m 



Carter et al. 2009 



the telescope and camera. The orbit-long slope probably 
arises from the combination of pointing/focus drifts (see 
Fig. [3]) with our imperfect flat-fielding of the detector. 
This effect of this orbital phase term could be equally 
well-achieved, for instance, by including a linear func- 
tion of the 0*^ order x and y positions ( see [Burke et al 
2010||Pont et al.|2007[jSwairi et al.|2008[ ). The visit-long 
slope is not mirrored in any of the measured geometri- 
cal properties of the star on the detector, and is more 
difficult to associate with a known physical cause. 

In order to determine the parameters C, V, B, R, D, 
and r, we fit Eq. [6] multiplied by a transit model to the 
last three orbits of each visit's uncorrected white light 
curve. The transit parameters are allowed to vary exactly 
as in §4.5.2[ including the use of the informative prior on 
the limb-darkening coefficients. The white light curves 
with the best model-ramp fit are shown in Fig. [51 and 
the properties of the residuals from this model are shown 
in Fig. [6| The transit parameters from this independent 
systematics correction method are consistent with those 
in Tab[3| We do not quite achieve the 280 ppm predicted 
scatter in the model-ramp light curves, and the residuals 
show slight evidence for correlated noise (Fig. [6| . More 
complicated instrumental correction models could almost 
certainly improve this, but we only present this simple 
model for heuristic purposes. In all sections except this 
one, we use the divide-oot-corrccted data exclusively 
for drawing scientific conclusions about GJ1214b. 

Fig. [8[ shows the inferred PDF's of the instrumental 
systematics parameters for all three HST visits, graphi- 
cally demonstrating the striking repeatability of the sys- 
tematics. As expected from the nearly identical cadence 
of illumination within each of the three visits, the ramp 
has the same R = 0.4% amplitude, r = 30 second 
timescale, and D = 20 second delay time across all ob- 
servations. The values of r and D are similar to the time 
for a single exposure, 25 seconds (including overhead). 
While the visit-long slope V is of an amplitude (fading 
by 0.06% over an entire visit) that could conceivably be 
consistent with stellar variability, the fact that it is iden- 
tical across all three visits argues strongly in favor of 
it being an instrumental systematic. B is the only pa- 
rameter that shows any evidence for variability between 
orbits; we would expect this to be the case if this term 
arises out of flat-fielding errors, since the 1*^* order spec- 
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Figure 8. The a posteriori distribution of ttie instrumental sys- 
tematics parameters from the analytic model, in each of the three 
visits. The MCMC results for single parameters {diagonal; his- 
tograms) and pairs of parameters (off-diagonal; contours encom- 
passing 68% and 95% of the distribution) are shown, marginalized 
over all other parameters (including c and d with priors, ti4, and 
-R*). V is measured in units of relative fiux/(3 X 96 minutes), B 
in relative flux/(96 minutes), R in relative flux, and both r and 
D in seconds. All visits are plotted on the same scale; for quanti- 
tative comparison, the median values and la uncertainties of each 
parameter are quoted along the diagonal. The systematics param- 
eters are remarkably repeatable from visit to visit; also, they are 
largely uncorrelated with Rp/Ri, [left column). 

trum falls on different pixels in the three visits. 

4.6. Spectroscopic Light Curve Fits 

We construct multiwavelength spectroscopic light 
curves by binning the extracted first order spectra into 
channels that are 5 pixels (AA = 23 nm) wide. We esti- 
mate the flux, flux uncertainty, and effective wavelength 
of each bin from the inverse- variance (estimated from the 
noise model) weighted average of each quantity over the 
binned pixels. For each of these binned spectroscopic 
light curves, we employ the divide-oot method to cor- 
rect for the instrument systematics. 

To measure the transmission spectrum of GJ1214b, we 
fit each of these 24 spectroscopic light curves from each 
of the three visits with a model in which Rp/ R^,, c, d, and 
s are allowed to vary. We hold the remaining parameters 
fixed so that a/R^ = 14.9749 and b = 0.27729, which 



are the valu es used by Bean et aL] ( 2010 ) , [Desert et al 



(201 la I, and [Croll et ar i 2Ull[ ). I'or limb-darkenmg pri- 
ors, we use the same sized Gaus sians o n the same linear 
combinations of c and d as in |4.5.1[ but center them 
on the PHOENIX-detcrmined best values for each spec- 
troscopic bin (see ^4.2). The correlation length of all 
parameters is < 10 m the MCMC chains. 

For most spectroscopic bins, the inferred value of s is 
within 1(7 of unity, indicating that the flux residuals show 
scatter commensurate with that predicted from photon 
noise (1400 to 1900 ppm across wavelengths). No evi- 
dence for correlated noise is seen in any of the bins, as 
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judged bythe same criterion as for the white light curves 
(see Fig^^[6). 

Fig. [9] snows the transmission spectra inferred from 
each 01 the three visits, as well as the divide-oot- 
corrected, spectrophotometric light curves from which 
they were derived. For the final transmission spectrum 
(shown as black points in Fig. |9]), we combine the three 
values of Rp/Ri,, and ctr^/r^ in each wavelength bin by 
averaging them over the visits with a weighting propor- 
tional to 1/ cr^ . Table |4j gives this average trans- 
mission spectrum, as well as the central values of the 
limb-darkening prior used in each bin. The wavelength 
grids in the three visits are offset slightly (by less than a 
pixel) from one another; in Table [4] we quote the average 
wavel ength for each bin. 



In i 



4.5.2| we found that GJ1214's 1% variability at 
WFC3 wavelengths causes AD = 0.014% or ARp/R^ = 
0.0006 variations in the absolute transit depth. The 
starspots causing this variability would have a similar 
effect on measurements of the transmission spectrum, 
but unless GJ1214's starspot spectrum is maliciously be- 
haved, the offsets should be broad-band and the influence 
on the wavelength-to-wavelength variations within the 
WFC3 transmission spectrum should be much smaller. 
Each visit's transmission spectrum is a differential mea- 
surement made with respect to the integrated stellar 
spectrum at each epoch; by averaging together three es- 
timates to produce our final transmission spectrum, we 
average over the time- variable influence of the starspots. 
Importantly, if GJ1214 is host to a large population of 
starspots that are symmetrically distributed around the 
star and do not appear contribute to the observed flux 
variability over the stellar rotation period, their effect 
on the transmission spectrum will not average out (se e 



Desert et al |2011b[ [Carter et al.|20Tl| [Berta et al.|2011[ ) 

It we fix the limb-darkening coetticients to the PHOENIX 
values instead of using the prior, the uncertainties on the 
Rp/Ri, measurements decrease by 20%. If we use only 
a single pair of LD coefficients (those for the white light 
curve) instead of those matched to the individual wave- 
length bins, the transmission spectrum changes by about 
1(7 on the individual bins, in the direction of showing 
stronger water features and being less consistent with an 
achromatic transit depth. These tests confirm that the 
presence of the broad H2O feature in the stellar spectrum 
(see Fig. p| makes it especially crucial that we employ 
the detailed, multiwavelength LD treatment. 

As a test to probe the infiuence of the divide-oot sys- 
tematics correction, we repeat this section's analysis us- 
ing the analytic model-ramp method to remove the in- 
strumental systematics; every point in the transmission 
spectrum changes by much less than Itr. We also exper- 
imented with combining the three visits' spectroscopic 
light curves and fitting for them jointly, instead of averag- 
ing together the transmission spectra inferred separately 
from each visit. We found the results to be practically 
identical to those quoted here. 

Because the transmission spectrum is conditional on 
the orbital parameters we held fixed (a/i?*, &), we under- 
estimate the uncertainty in the absolute values of Rp/ Ri, ; 
the quoted (Jr^/r^^ are intended for relative comparisons 
only. Judging by the Rp/Ri, uncertainty in the uncon- 
strained white light curve fit (Table [2]), varying a/R^, 



could cause the ensemble of Rp/R^, measurements in Ta- 
ble[4|to move up or down in tandem with a systematic un- 
certainty that is comparable to the statistical uncertainty 
on each. This is in addition to the A Rp/Rj, = 0.0006 off - 
sets expected from stellar variability ( Berta et al.|2011 ). 



4.7. Searching for Transiting Moons 

Finally, we search for evidence of transiting satellite 
compansions to GJ1214b in our summed WFC3 light 
curves. The light curve morphology of transiting exo- 
moons can be complicated, but they could generally ap- 
pear in our data as shallow transit-shaped dimmings or 
brightenings offset from the planet's transit light curve 
(see Kipping 2011, for a detailed discussion). While the 
presence ol aTmoon could also be detected i n tempo- 
ral va riations of the planetary transit duration ( Kipping 
20091, we only poorly constrain GJ1214b's transit dura- 



tion m individual visits due to incomplete coverage. 

Based on the Hill stability criterion, we would not ex- 
pect moons to survive farther than 8 planetary radii away 
from GJ1214b so their transits should not be offset from 
GJ1214b's by more than 25 minutes, less than the dura- 
tion of an HST Earth occultation. We search only the 
data in the in-transit visit, using the divide-oot method 
to correct for the systematics. Owing to the long buffer 
download gaps in our light curves (see q2]), the most 
likely indication of a transiting moon in theWFCS light 
curve would be an offset in flux from one 12-exposure 
batch to another. Given the 376 ppm per-exposure scat- 
ter in the divide-oot corrected light curve, we would 
have expected to be able to identify transits of 0.4 
R® (Ganymede-sized) moons at 3(7 confidence. We see 
no strong evidence for such an offset. Also, we note that 
starspot occultations could easily mimic the light curve 
of a transiting exomoon in the time coverage we achieve 
with WFC3, and such occ ultations are known to occur in 



the G J1214b system (see [j3erta et al.[[2011[ [Carter et al. 
2011||Kundurthy et al.[[20TTf : 



Due to the many possible configurations of transiting 
exomoons and the large gaps in our WFC3 light curve, 
our non-detection of moons does not by itself place strict 
limits on the presence of exo-moons around GJ1214b. 

5. DISCUSSION 

The average transmission spectrum of GJ1214b from 
our three HST visits is shown in Fig. [9[ To the precision 
afforded by the data, this transmission spectrum is fiat; 
a simple weighted mean of the spectrum is a good fit, 
with = 20.4 for 23 degrees of freedom. 

5.1. Implications for Atmospheric Compositions 

We compare the WFC3 transmission spectrum to a 
suite of cloud-free theoretical atmosph ere models for 
GJ1214b. The models were calculated in IMiller-Ricci fcl 
Fortney[ ( 2010|, and we refer the reader to that paper for 
their details, 'i'o compare them to our transmission spec- 
trum, we bin these high-resolution (i? — 1000) models to 
the effective wavelengths of the 5-pixel WFC3 spectro- 
scopic channels (i? = 50 — 70) by integrating over each 
bin. Generally, to account for the possible suppression 
of transmission spectrum features caused by the overlap 
of shared planetary and stellar absorption lines, this bin- 
ning should be weighted by the photons detected from 



WFC3 Observations of GJ1214b 



13 



-40 
-20 

20 
40 
-40 

-20 

20 
40 



1 .45 
S 1.40 

1 35 
^ 1.30 

1.25 





after 
divide-oot - 
correction 




Visit 1 
I Visit 2 

Visit 3 
combined 



1.1 



1.2 



1.4 1.5 
Wavelength (/U,m) 



1.6 



1 .7 



Figure 9. Top panels: Spectroscopic transit light curves for GJ1214b, before and after the divide-oot correction, rotated and offset for 
clarity. Bottom panel: The combined transmission spectrum of GJ1214b {black circles with error bars), along with the spectra measured 
for each visit (colorful circles). Each light curve in the top panel is aligned to its respective wavelength bin in the panel below. Colors 
denote HST visit throughout. 



the system at very high resolution, but this added com- 
plexity is not justified for our dataset. The normalization 
of the model spectra is uncertain (i.e. the planet's true 
Rp), so we allow a multiplicative factor in Rp/R^, to be 
applied to each (giving 24-1=23 degrees of freedom for 
all models). Varying the bin size between 2 and 50 pixels 
wide does not significantly change any of the results we 
quote in this section. 

A solar composition atmosphere in thermochemical 
equilibrium is a terrible fit to the WFC3 spectrum; it 
has a x^=126.2 (see Fig. ^ and is formally ruled out 
at 8.2(7 confidence. Likewise, the same atmosphere but 
enhanced 50 x in elements heavier than helium, a quali- 
tative approximation to the metal enhancement in the 
Solar System ice giants (enhanced 30 — 50x in C/H; 



Gaut ier et al 
rulei 



1995 



Encrenazlf2005 

IF 



Guill ot fc Gautier 
BotE 



2009P , is ruled out at 7.5a (x^ = 113.2). Both models 
assume equilibrium molecular abundances and the ab- 
sence of high-altitude clouds; if GJ1214b has an H2-rich 
atmosphere, at least one of these assumptions would have 
to be broken. 

Suggesting, along these lines, that photoch emistry 
might deplete GJ1214b's atmosph ere of methane, ' Desert 
et^ ([20TTal), |Croll et _al| (j20lT|, and [Crossficld et al. 
I^UTT) have noted their observations to be consistent 
with a solar composition model in which CH4 has been 
artificially removed. With the WFC3 spectrum alone, 
we can rule out such an H2-rich, CH4-free atmosphere 



Kempton et aT 



at 6.1a (Fig. 10 ) . Thi s is consistent with Miller- Ricci 



( 201l"] )'s theoretical finding that such 
thorough methane depletion cannot be achieved through 
photochemical processes, even when making extreme as- 



sumptions for the photoionizing UV flux from the star. 

P revious spectroscopic m easurements in the red opti- 
cal ( Bean et al.|2010 2011 ) could only be reconciled with 
a H2-rich atmosphere it such an atmosphere were to host 
a su bstantial cloud layer at an altitud e above 200 mbar 
(see'Mille r-Ricci Kempton et al.|2011 1. How far the flat- 
tening influence of such a cloud layer would extend be- 
yond 1 /im to WFC3 wavelengths would depend on both 
the concentration and size distribution of the scattering 
particles. As such, we explore possible cloud scenarios 
consistent with the WFC3 spectrum in an ad hoc fash- 
ion, using a solar composition atmosphere and arbitrarily 
cutting off transmission below various pressures to emu- 
late optically thick cloud decks at different altitudes in 
the atmosphere. Fig. [Tl] summarizes the results. A cloud 
deck at 100 mbar, which would be sufficient to flatten the 
red optical spectrum, is ruled out at 5.7a {x^ = 82.8). 
Due to higher opacities between 1.1 and 1.7 /im, WFC3 
probes higher altitudes in the atmosphere than the red 
optical, requiring clouds closer to 10 mbar to match the 
data (x^ = 23.4). Note, with the term "clouds" we refer 
to all types of particles that cause broad-band extinction, 
whether they scatter or absorb, and whether they were 
formed through near-equilibrium condensation (such as 
Earth's water clouds) or through upper atmosphere pho- 
t ochemistry (suc h as Titan's haze) 



Kcr 



Fortneyl (|2005[) and [Miller- Ricci Kempton et al. ( |2011 | 
identified 



and ZnS as condensates that would be 
likely to form in GJ1214b's atmosphere, but found they 
would condense deeper in the atmosphere (200-500 mbar) 
than required by the WFC3 spectrum and would proba- 
bly not be optically thick. While winds may be able to 
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Figure 10. The WFC3 transmission spectrum of GJ1214b {black circles with error bars) compared to tlieoretical models {colorful lines) 
witli a variety of compositions. The high resolution models are shown here smoothed for clarity, but were binned over each measured 
spectroscopic bin for the comparisons. The amplitude of features in the model transmission spectra increases as the mean molecular 
weight decreases between a 100% water atmosphere {fi = 18) and a solar composition atmosphere {fi = 2.36). 
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Figure 11. The WFC3 transmission spectrum of GJ1214b {black 
circles with error bars) compared to a model solar composition 
atmosphere that has thick clouds located at altitudes of 100 mbar 
{pink lines) and 10 mbar red lines). We treat the hypothetical 
clouds in an ad hoc fashion, simply cutting off transmission through 
that atmosphere below the denoted pressures. 



loft such clouds to higher altitudes, it is not clear that 
the abundance of these species alone would be sufficient 
to blanket the entire limb of the planet with optically 
thick clouds. The condensation and complicated evo- 
lution of clouds has been studi ed within the context of 
cool sta rs an d hot Jupiters (e.g|Lodders & Fegley 200^ 
[Helling et al. 2008 ) , but further study uito the theoretical 
landscape for equilibrium clouds on planets in GJ1214b's 
gravity and temperature regime is certainly warranted. 
The scattering may also be due to a high altitude haze 
formed as by-products of high-alti tude photochemistry; 



Miller- Ricci Kempton et al. (2011) found the conditions 
on GJ 1214b to allow for tfie forrnation of complex hy- 
drocarbon clouds through methane photolysis. 

However such clouds might form, they would either 
need to be optically thick up to a well-defined altitude 
or consist of a substantial distribution of particles act- 
ing in the Mie regime, i.e. with sizes approaching 1 
fim. Neither the VLT spectra nor our observations give 
any definitive indications of the smooth falloff in transit 
depth toward longer wavelengths that would be expected 
from Rayleigh scattering by molecules or small particles. 
This is unlike the case of the hot Jupiter HD189733b, 
where the uniform decrease in transit depth from 0.3 to 
1 liui dP ont et al.l2008; Leca velier Pes Etangs et al.|2008"' 



Sing et al. 2011) and perhaps to as tar as 3.(j /nm (see |!Siiig 
erarf2009, ^Desert et al.^2009| has been convincingly at 



tributed to a small particle haze . 

As an alternative, the transmission spectrum of 
GJ1214b could be flat simply because the atmosphere has 



a large mean molecular weight. We test this possibility 
with H2 atmospheres that contain increasing fractions of 
H2O. This is a toy model, but including molecules other 
than H2 or H2O in the atmosphere would serve prin- 
cipally to increase /i without substantially altering the 
opacity between 1.1 and 1.7 fim, so the limits we place 
on ^ are robust. We find that an atmosphere with a 
10% water by number (50% by mass) is disfavored by 
the WFC3 spectrum at 3.1tT (x^ = 47.8), as shown in 
Fig. |9] All fractions of water above 20% (70% by mass 
are good fits to the data (x^ < 25.5). The 10% wa- 
ter atmosphere would have a minimum mean molecular 
weight of /X = 3.6, which we take as a lower limit on the 
atmosphere's mean molecular weight. 

For the sake of placing the WFC3 transmission spec- 
trum in the context of other observations of GJ1214b, we 
also display it alongside the published tran smission spec 



tra from the VLT (Bean et al.|2010 



et all 12011 1), Magel 



( |Uesert et aL||2011a[ ) in Fig 



cause individual se 



an 



(IBean et al.lT 
FiK.[T2l SteE 



2011 1 CFHT ( [Croll 
|20li|), and Spitzer 



Stellar variability could 
of observations to move up and 



down on this plot by as much as AD = 0.0 14% for mea- 
surements in the near- IR ( [Berta et al.|2011[ ); we indicate 
this range of potential offsets by an arrow at the right of 
the plot. We display the measurements in Fig. 12 with 
no relative offsets applied and note that their general 
agreement is consistent with the predicted small influ- 
ence of stellar variability. Depending on the temperature 
contrast of the spots, however, the variability could be 
larger by a factor of 2 — 3x in the optical, and we cau- 
tion the reader to consider this systematic uncertainty 
when comparing depths between individual studies. For 
instance, the slight apparent rise in Rp/Ri, toward 0.6 
/xm, that would potentially be consistent with Rayleigh 
scattering in a low-/i atmosphere, could also be easily ex- 
plained through the poorly constrained behavi or of the 
star in the optical. 



Indeed, 



Bean et al. (2011) found a 



significant offset between datasets that overlap in wave 
length (near 0.8 /im) but were taken in different years, 
suggesting variability plays a non-negligible role at these 
wavelengths. 

Finally, we note that any model with /i > 4, such as 
one with a > 50% mass fraction of water, would be con - 

(pOlO ) 



sistent with the measurements from Bean et al 



Desert et alT] ( [2011a[ ), [Crossfield et all ( |201ip , [Bean et al 
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Table 4 

GJ1214b's Transmission Spectrum from WFC3/G141 







Rp/R* 
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1.123 





11641 ± 0.00102 


-0.372 


1.068 


1.146 





11707 ± 0.00099 


-0.397 


1.088 


1.169 





11526 ±0.00098 


-0.404 


1.089 


1.192 





11589 ±0.00093 


-0.410 


1.090 


1.215 





11537 ± 0.00091 


-0.406 


1.075 


1.239 





11574 ± 0.00090 


-0.403 


1.063 


1.262 





11662 ±0.00088 


-0.407 


1.064 


1.285 





11565 ±0.00088 


-0.403 


1.045 


1.308 





11674 ± 0.00085 


-0.411 


1.042 


1.331 





11595 ±0.00087 


-0.390 


1.068 


1.355 





11705 ± 0.00089 


-0.368 


1.101 


1.378 





11664 ± 0.00088 


-0.371 


1.110 


1.401 





11778 ± 0.00088 


-0.338 


1.075 


1.425 





11693 ±0.00091 


-0.295 


1.029 


1.448 





11772 ±0.00090 


-0.319 


1.056 


1.471 





11663 ± 0.00092 


-0.322 


1.061 


1.496 





11509 ±0.00100 


-0.345 


1.084 


1.517 





11635 ±0.00104 


-0.305 


1.013 


1.541 





11626 ±0.00091 


-0.330 


1.042 


1.564 





11681 ± 0.00091 


-0.351 


1.059 


1.587 
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( |2011 ), and WFC3. The only observation it could not 
explain wou ld be the deep iiTs-measurement from |Croll| 



et al.| |2011|. Of the theoretical models we tested, we 
could find none that matched all the available measure- 
ments. We are uncertain of how to interpret this appar- 
ent incompatibility but hopeful that future observational 
and theoretical studies of the GJ1214b system may clar- 
ify the issue. In the meantime, we adopt an atmosphere 
with at least 50% water by mass as the most plausible 
model to explain the WFC3 observations. 

5.2. Implications for GJ1214b's Internal Structure 

If GJ1214b is not shrouded in achromatically optically 
thick high-altitude clouds, the WFC3 transmission spec- 
trum disfavors any proposed bulk composition for the 
planet that relies on a substantial, unenriched, hydrogen 
envelope to explain the planet's large radius. Both the 
ice-rock core with nebular H/He envelope and pure rock 
core with outgas s ed H2 envelope scenarios explored by 
Rogers &: Seager (2010) would fall into this category, re- 
quiring additional mgredients to match the observations. 
In contrast, their model that achieves GJ1214b's large 
radius mostly from a large water-rich core, would agree 
with our observations. 

Perhaps most compellingly, a high /i sce nario would 
be consisten t with composition proposed by [Nettelmann 
et al. (2011), who found that GJ1214b's radius could be 
explained by a bulk composition consisting of an ice-rock 
core surrounded by a II/IIe/Il2 envelope that has a wa- 
ter mass fraction of 50-85%. Such a composition would 
be intermediate between the H/He- and H20-env elope 



hmiting cases proposed by Rogers fc Seager] (2010). The 
H/He/H2 envelope might arise it GJ12l4b had origi- 
nally accreted a substantial mass of hydrogen and helium 
from the primordial nebula but then was depleted of its 
lightest molecules through atmospheric escape. 

5.3. Prospects for GJ1214b 



Futu re observations with the James Webb Space Tele - 
scope ( [Deming et al]|2009[ p^altenegger fc Traub||2009[ ), 
one of the immense next gener ation ground- b ased tele- 
scopes (GMT, TMT, ELT; see [Ehrenreich et al. 2006), 
or possibly even a dedicated campaign with current fa- 
cilities, could detect the 0.01% transmission spectrum 
features of a 100% water atmosphere on GJ1214b, and 
potentially distinguish between clear H2-poor and cloudy 
H2 -rich atmosphe res. Along another front, simulations 



by Menou 
thermal p. 



(2011) show that observations of GJ1214b's 



lase curve, such as those for HD189733b by 



Knutson et al. (2007a), would probe the ratio of radiative 
to advective tiinescaies in GJ1214b's outer envelope and 
provide an independent constraint on the atmospheric 
composition. Detecting the thermal emission from this 
500K exoplanet is currently very difficult, and will likely 
have to wait until the launch of JWST. 

In the meantime, we advocate further study of the 
GJ1214 system in gener al. Confirming and ref ining the 
parallax for the system ( van Altena et al?]|1995p will im- 
prove our knowledge of the stellar mass, and in turn, the 
planet's mass and radius. Likewise, further radial veloc- 
ity ob servations wil l empi rically constrain the hypothe- 
sis by Carter et al. (2011) that a significantly non-zero 
orbital eccentricity could be biasing GJ1214b's inferred 
density. 

6. CONCLUSIONS 

In this work, we made new measurements of the 
GJ1214b's transmission spectrum using HST/WFC3. 
Reaching a precision of tr^/r^ = 0.0009 in 24 simultane- 
ously measured wavelength bins, we found the transmis- 
sion spectrum to be completely fiat between 1.1 and 1.7 
/im. We saw no evidence for the strong H2O absorption 
features expected from a range of H2-rich model atmo- 
spheres. 

Given the lack of a known source for clouds or hazes 
that could create a truly achromatic transit depth across 
all wavelengths, we interpret this fiat WFC3 transmission 
spectrum to be best explained by an atmosphere with a 
high mean molecular weight. Based on our observations, 
this atmosphere would likely consist of more than 50% 
water by mass or a mean molecular weight of /.t > 4. Such 
an atmosphere would be consistent with observa t ions o f 
GJ1214b's transmiss i on spectrum by Bean et aT] (| 2010 ) 



Desert et all ( |2011a[ ), fCrossfield et alj (|201l[ ),"an(l 'Bean| 
et al. ( 2011[ ) although it would be dimcult to reconcile 



wIHrthosTby lCroll et ah] ([20ll| 

Such a constraint on GJ12l4b's upper atmosphere 
serves as a boundary condition for models of bulk com- 
position and structure of the rest of the planet. It sug- 
gests GJ1214b contains a substantial fraction of water 
throughout the interior of the planet in order to obviate 
the need for a completely H/He- or H2-dominatcd enve- 
lope to explain the planet's large radius. A high bulk 
volatile content would point to GJ1214b forming beyond 
the snow line and migrating inward, although any such 
statements about GJ1214b's past are subject to large 
uncertainties in th e atmospheric mass loss history (see 
[Rogers et al.|2011| ). 

Finally, this paper is the first published study using 
WFC3 for observing a transiting exoplanet. Aside from 
several instrumental systematics that were straightfor- 
ward to correct and did not require a detailed instrumen- 
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Figure 12. GJ1214b's transmission spec trum from WFC3 in th e context of obser vations from the VLT (0.6-1 lira; IBean et al.|2010 |[201l1l, 
CFHT (1.25 + 2.15 /^m; [Croll et al.|2011[ l, Magellan (2.0-2.3 ^m |Bean et al.l2011| , and Spitzer (3.6 -I- 4.5 Ajm;'Deser! et al.l201i^ While 
they do not measure an absolute transit depth, observations from INiKBJr'JDC on Keck (2.1-2.3 /im Crossficld et al. 2011) disfavor models 
they tested that had amplitudes larger than 0.05% in their wavelength range; we represent these constraints with the dashed rectangular 
boxes. Two extr emes of the models explored in this paper are shown, normalized to the MEarth-measured transit depth (see Millcr-Ricci| 

jcn 



\&L Fortney|20ldll . It is important to note that stellar variability could cause individual data sets to shift up or down on this plot as mucB 
as l\U = J.U1470 in the near-IR or 2 — 3x more in the optical, depending on the stellar spot spectrum. 



tal model, the camera delivered nearly photon-limited 
performance both in individual spectrophotometric light 
curves and in summed white light curves. We are con- 
fident that WFC3 will serve as a valuable tool for exo- 
planet atmospheric characterization in the years to come. 
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